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We study electron transport through a domain wall in a ferromagnetic nanowire subject to spin- 
dependent scattering. A scattering matrix formalism is developed to address both coherent and 
incoherent transport properties. The coherent case corresponds to elastic scattering by static de- 
fects, which is dominant at low temperatures, while the incoherent case provides a phenomenological 
description of the inelastic scattering present in real physical systems at room temperature. It is 
found that disorder scattering increases the amount of spin-mixing of transmitted electrons, reduc- 
ing the adiabaticity. This leads, in the incoherent case, to a reduction of conductance through the 
domain wall as compared to a uniformly magnetized region which is similar to the giant magnetore- 
sistance effect. In the coherent case, a reduction of weak localization, together with a suppression 
of spin-reversing scattering amplitudes, leads to an enhancement of conductance due to the domain 
wall in the regime of strong disorder. The total effect of a domain wall on the conductance of a 
nanowire is studied by incorporating the disordered regions on either side of the wall. It is found 
that spin-dependent scattering in these regions increases the domain wall magnetoconductance as 
compared to the effect found by considering only the scattering inside the wall. This increase is 
most dramatic in the narrow wall limit, but remains significant for wide walls. 



I. INTRODUCTION 

The interplay between magnetic structure and elec- 
trical resistance in mesoscopic ferromagnetic systems 
is interesting from both technological and fundamental 
points of view. Magnetoresistance properties of domain 
walls represent a particularly striking example of this 
rich physical problem, and have been intensively stud- 
ied in recent years. Experiments on ferromagnetic thin 
films originally found that domain walls contribute to 
an enhancement of conductance, 1,2 although this was 
later understood to be due to anisotropic magnetoresis- 
tance (AMR) effects, 3 Further experiments on cylindrical 
Co nanowires^ and thin polycrystalline Co films^ have 
yielded a reduction in conductance, distinct from the 
positive AMR contribution, which has been attributed 
to scattering from magnetic domain walls. 

From the theoretical side two main approaches exist 
to address this problem. On one hand, first-principle 
calculations take into account realistic band structures 
to study equilibrium properties. Assuming that the re- 
sulting one-body wavefunctions are a good description of 
the many-body wavefunctions, a calculational scheme to 
study electronic transport can be developed. According 
to the details of the model, large domain wall magne- 
toresistances may be found&i On the other hand, phe- 
nomenological models including only the essential fea- 
tures of the band structure are extremely useful as they 
give insight into important physical mechanisms. The 
most widely used approach is the sd model, originally 
due to Mottj£ which separates the conduction (s) elec- 
trons from those responsible for the magnetic structure 
(rf). The exchange interaction between s and d electrons 
is incorporated in the spin-splitting A of the conduction 
band. An important advantage of this model is that it 



can be readily extended to include impurity scattering 
through a disorder potential, making it particularly suit- 
able for the topic of the present work. 

The simplest approach is to treat transport through 
the domain wall and the leads as ballistic, where the only 
scattering is due to the rotating spin-dependent poten- 
tial in the wall. The resulting magnetoconductance is 
very small for typical domain wall widthsj2iifi although 
it becomes significant in the special case of domain walls 
trapped in magnetic nanocontacts>ii*i2ii^ In general, the 
most important effect of a domain wall on ballistic elec- 
tron transport is a mixing of the up and down spin chan- 
nels, which arises from the inability of electrons to follow 
adiabatically the local magnetization direction, referred 
to as "mistracking" >Mii^ When spin-dependent scatter- 
ing in the regions adjacent to the domain wall is taken 
into account^ this mistracking of spin in a ballistic wall 
leads to a significant magnetoconductance analogous to 
the giant magnetoresistance effect^ 

However, treating the domain wall as a completely bal- 
listic (disorder-free) system is not realistic for experimen- 
tally relevant systems. For instance, the cobalt nanowires 
of Ref. have an estimated elastic mean free path of 
^7nm, while the wall width in cobalt is ~15nm. Since 
these two characteristic length scales are of the same or- 
der, the transport through the domain wall cannot be 
described as either ballistic or diffusive. For other ma- 
terials, such as iron or nickel, the wall widths are larger 
and, depending on the amount of disorder, the diffusive 
regime may be reached. 

A number of works have focused on the role of disor- 
der scattering inside the domain wall. Viret et ali^ used 
the ballistic mistracking of spin mentioned above to de- 
velop an intuitive picture based on a weighted average of 
up and down resistances. The estimated relative magne- 
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toresistance decreases with the width of the domain wall 
and agrees with measurements on domain wall arrays in 
thin films. Similar results were found using models based 
on the Boltzmann equation^ and the Kubo formula^* 
Phase coherence effects have also been studied^LSi and 
in the case of spin-independent disorder a negative con- 
tribution to domain wall magnetoconductance has been 
predicted. 20 ' 22 

In this work we develop a new model for transport 
through a disordered domain wall based on combin- 
ing scattering matrices for individual impurity scatter- 
ers, which improves on existing treatments in several 
key aspects. Firstly, the model is non-perturbative in 
wall width and disorder strength, which allows us to 
study walls of arbitrary width and consider both ballistic 
and diffusive transport regimes. It is found that impu- 
rity scattering inside the domain wall causes an increase 
in transmission and reflection with spin-mistracking, or 
equivalently a reduction in the adiabaticity of spin trans- 
port through the wall. Secondly, we can treat both 
phase-coherent and incoherent transport regimes within 
the same model. This permits a quantitative determi- 
nation of the contribution of phase coherence effects to 
domain wall magnetoconductance. In this way we find 
that for incoherent transport a domain wall gives rise 
to a positive magnetoconductance effect {i.e. a reduction 
of conductance), which depends on the relative impurity 
scattering strength for up and down electrons and scales 
linearly with the number of conductance channels. In 
constrast, for coherent transport in the case of strong dis- 
order, the domain wall magnetoconductance is negative 
and does not depend on the relative up/down scatter- 
ing strength or the system size. Finally, spin-dependent 
scattering in the regions adjacent to the domain wall can 
be incorporated directly into the model, allowing us to 
calculate the total magnetoconductance effect of a do- 
main wall in a nanowire using an approach similar to the 
circuit model developed in Ref. It is found that scat- 
tering in the uniformly magnetized regions on either side 
of the wall causes an enhancement of the magnetocon- 
ductance effect which is largest for narrow walls but is 
also significant in the wide wall limit. 

The layout of this paper is as follows. In Section II we 
introduce our physical model and describe the numerical 
method which we use to calculate conductance through 
a disordered region. In Section III we study the intrinsic 
transport properties of a disordered domain wall in both 
coherent and incoherent regimes. In Section IV we incor- 
porate spin-dependent scattering in the regions adjacent 
to the wall. Finally, in Section V we discuss the exper- 
imental relevance of our findings, presenting our conclu- 
sions and outlook. 



II. PHYSICAL MODEL 

We consider a quasi one-dimensional (quasi-lD) wire, 
with longitudinal axis lying along the z axis and a sin- 



gle transverse dimension of length L y (see Figure^ . The 
position r denotes the two-dimensional vector (y, z). The 
extension to a two-dimensional cross section is straight- 
forward, but complicates the notation and numerical 
calculations without adding new physics. Within the 
sd model, the conduction electrons in the wire are de- 
scribed by a free-electron one-body hamiltonian with 
spin-dependent potential 

H = -^ 2 + jM-a + V(f). (1) 

Here /(f) is a unit vector in the direction of the effective 
field representing the magnetic moment due to the d elec- 
trons, a is the Pauli spin operator of the s electrons, and 
A is the spin-split energy gap between up and down s 
electrons. The potential V(f) represents impurities and 
leads to scattering. It is discussed in detail in Section 

urn 

For a wire uniformly magnetized in the z direction we 
have /(f) = (0, 0, —1) for all r. When the wire contains a 
domain wall separating regions of opposite magnetization 
along the z axis, we have /(f) = (0, 0, —1) for z < and 
/(f) = (0, 0, 1) for z > A, with the length A defining the 
wall region < z < A. 

For simplicity we assume a square well potential with 
infinite walls at y = 0, L y for the transverse confinement. 
In the disorder-free regions of constant magnetization at 
either end of the system (i.e. the leads), electrons then 
occupy well-defined modes (channels): 

< 2 > 

where n is a positive integer. 

The eigenstates in the leads (which constitute the 
asymptotic states in a scattering approach to the do- 
main wall) are characterized by a longitudinal wavevec- 
tor k, transverse mode number n, and spin eigenvalue 
a = ± (representing spin states which are, respectively, 
anti-parallel and parallel to the local direction of /(f)). 
The corresponding dispersion relation is 

E _bW /mr\ 2 aA 

The relevant states for transport are those with energy E 
at the Fermi energy E-p , and corresponding wavevectors 

Here fc^p — \/k F + ak\j2 is the spin-dependent Fermi 
wavevector, defined in terms of k-p = \J 2mEp / h 2 and 
k& = y^2mA/h 2 . The number of propagating modes in 
each lead and for each spin sub-band, N a , is given by the 
largest value of n for which k an is real (N a ~ h^pLy/ir). 
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The wavefunctions are two-component spinors: 

V4(r) = x [-^e ±ik ^Uy) W) , (5) 

V mv<j n 

where ^ denotes propagation to the right /left, and the 
spinor basis states |±) represent spin eigenstates parallel 
or anti-parallel to the local value of f(r). The states are 
normalized to unit flux, where the velocity factor for a 
plane-wave state is v an = hk an /m. 

Defining the spin states \a) with respect to the local 
magnetization direction allows us to treat both of the 
leads (as well as the uniform case) within the same nota- 
tion. Inside the wall (0 < z < A), \a) depend on r, and 
hence the relevant electron eigenstates are more compli- 
cated than the states in Eq. (J5J. They are described in 
Section ITTBl 

The parameters Ep and A characterize the structure 
of the parabolic s band in the sd model, and must be cho- 
sen to approximate the spin polarization obtained from 
band structure calculations. This approximation intro- 
duces considerable uncertainty. In this work we find 
that the dependence on Ep and A, as well as the wall 
width A, can be encapsulated in a dimensionless "effec- 
tive" wall width pp (Eq. I^bj)). Throughout this work, 
we fix Ep — lOeV and A = O.leV, which correspond to 
reasonable parameter estimates for cobalt^ By varying 
A, we then obtain results which correspond to different 
effective domain wall widths and apply to arbitrary ma- 
terials. 



A. Scattering matrix approach to conductance 

Since we assume non-interacting electrons, the con- 
ductance g (in units of e 2 /h) may be described by the 
Landauer-Buttiker formula^: 

( N„, N„ ] 

9= E \ EEww , (6) 

er',<7=± I n' = ln=l I 

where T (T i n i :an (Ep) is the probability of an electron at 
energy Ep to be transmitted from the state with spin 
and transverse mode (a, n) in the left lead to the one 
with (cr',n') in the right lead. 

To calculate the transmission amplitudes and probabil- 
ities, we use an approach based on combining the scat- 
tering matrices of all the scatterers in the system, which 
was developed to treat transport in spin-independent dis- 
ordered systems by Cahay et alm^ 

Our composition of scattering matrices neglects the ef- 
fect of evanescent modes, which have imaginary wavevec- 
tor and are assumed to decay between successive scatter- 
ing events. This approximation has been widely used in 
the so-called local approach for the random matrix de- 
scription of quasi- ID wires iSSiSL2& The impressive success 
of random matrix theory in explaining the universal fea- 
tures of quantum transport, as well as its agreement with 
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FIG. 1: Schematic illustration of the scattering processes oc- 
curring in our model of a disordered wire with (a) uniform 
magnetization and (b) a domain wall. Above and below the 
wire are arrows indicating the magnetization direction, f(r), 
as well as the scattering matrices describing the transport 
through the wire in each case. The impurities are represented 
by point scatterers randomly located at f a (a = 1, . . . ,N[), 
with corresponding scattering matrices s a (uniform) and s Q 
(domain wall) . The matrices p a and p~ a represent the ballistic 
propagation between successive scatterers, while s( L,R ) repre- 
sent the scattering from the domain wall interfaces (case (b) 
only). These matrices can be combined in a coherent or inco- 
herent way to give the total scattering matrix for the system 
(Eqs. 1121 and 1131 respectively). The incident and scattered 
fluxes, 

ttjn i kin and a out , bout, are also shown. 



numerical simulations and microscopic theories, relies on 
the hypothesis of a quasi-lD geometry (L z 3> L y ) and 
weak scattering (I ^> SL ^> kp~ , where 5L is the size of 
the scattering blocks). In our case these two hypotheses 
are valid, since the mean distance between scatterers is 
much larger than the Fermi wavelength. 

For the problem of spin-dependent transport our 
method offers several important advantages over the 
spin-dependent extension of the recursive Green's func- 
tion (RGF) technique. The latter calculates the phase- 
coherent conductance with a tight-binding approxima- 
tion for Eq. Q in which the rotating magnetization and 
the impurity potential are incorporated through spin- 
dependent on-site potential energies. 15,29 In particular, 
the scattering matrix approach allows the conductance to 
be calculated either coherently or incoherently, which is 
useful for identifying and understanding phenomena aris- 
ing from phase coherence. Additionally, because scatter- 
ing matrices are expressed in terms of local basis states, 
we are able to avoid the coarse-grained discretization in- 
herent in the tight-binding description of the domain wall 
magnetization. Finally, because only propagating states 
are included in the scattering matrices, the computa- 
tional effort required for a system with a given number 
of conducting channels is significantly lower than with 
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the RGF technique. In the case of coherent transport 
through a disordered domain wall, we have verified that 
results using the scattering matrix approach are in quan- 
titative agreement with those obtained using the RGF 
technique^ 

The scattering matrix of a system relates incoming 
(flin, kin) and outgoing (a ou t, b out ) flux amplitudes from 
the left and right (Figure through 



frout 
flout 



t r' 

r a 



(7) 



Here t, r (t',r') are transmission and reflection matrices 
for states incident from the left (right). The elements 
of these sub-matrices are the scattering amplitudes be- 
tween individual modes (cr, n) and (cr', n'), which we write 
generically as £a'n';<rn> where £ stands for one of t, r, t', 
or r' . The dimension of these sub-matrices of s is deter- 
mined by the number of propagating modes of each spin 
on either side of the scatterer. 

Throughout this work we will denote amplitudes and 
their matrices by lower-case letters and the corresponding 
probabilities (given by the absolute square of the ampli- 
tudes) by upper-case letters. We therefore use s to denote 
scattering matrices containing amplitudes and S to de- 
note the corresponding matrix of probabilities. This con- 
vention differs from the standard notation in which 5* de- 
notes the scattering matrix of amplitudes, but we adopt 
it here for notational consistency. The sub-matrices of S 
are the transmission and reflection probability matrices 
T, R, T', and R' . We refer to these matrices generically 
as S, and the corresponding individual probabilities as 
^a'n'-an = \£,a'n'-an \ 2 ■ We will also make use of the fol- 
lowing notation for sums of transmission and reflection 
probabilities summed over transverse modes and spin: 



n'n 



■ j G'n' van > 



(8a) 



(8b) 



Here S CT > CT represents the total scattering probability from 
spin sub-band cr to cr', and S CT the total scattering proba- 
bility (into both spin sub-bands) for states incident with 
spin cr. 

For two scatterers described by the individual scat- 
tering matrices s\ and s 2 , the resultant scattering ma- 
trix obtained by combining them in series is written as 
s\2 = Si * s 2 , where * stands for the composition law 

tu = fa • (1 - r[ ■ r 2 )~ l • t u (9a) 

ru = n +t[ T 2 (l-ri T2)- 1 (9b) 

*i 2 = ti(l + r 2 (l - r[ ■ r 2 )- x ■ r[) ■ t' 2 , (9c) 

r' u = r' 2 + 1 2 ■ (1 - ri • r^ 1 ■ r[ ■ t' 2 . (9d) 

Here 1 denotes the identity matrix with the same number 
of rows as r[ and the same number of columns as r 2 . 
The application of this composition law to include many 
scatterers forms the basis of the approach used in this 
work. 



1. Coherent and incoherent transport through a disordered 
region 

In this work we are interested in the conductance 
through a disordered region with either a domain wall 
or uniform magnetization. As we shall discuss in Section 
III CI the effect of disorder is represented in our model by 
the potential V(r), which is comprised of N[ delta func- 
tion scatterers at randomly distributed positions f a . The 
indices are ordered so that z a < z a+ \ for a = 1, . . . , JVj. 

Within the scattering matrix approach, the phase- 
coherent transmission from z = to A is calculated by 
combining the scattering matrices of each delta function 
scatterer with those for ballistic propagation between the 
scatterers. We denote these matrices respectively by s a , 
p a for uniform magnetization, and s a , p a for a domain 
wall. In the case of a domain wall, we also require scat- 
tering matrices s( L,R ' to describe the scattering at the 
interfaces between the domain wall and uniform regions. 
These matrices, and the basis states inside the domain 
wall, are discussed in detail in the next sub-section and 
in Appendix IA1 

The matrices p a and p a contain the phase shifts 
acquired by electrons propagating ballistically between 
scatterers at r a _i and f a . We write the longitudi- 
nal propagation distances as 8z a = z a — z a -i, where 
z±, . . . zn { are the longitudinal components of the impu- 
rity positions and we define zq — 0, z/v i+ i = A. For the 
uniform case we have 



Pa 



where 



cr' ri :crn 



(f a 

tp a 



fia'aSn'n^X'p {lk an 8z 



(10) 



(11) 



The matrices p a for the domain wall case are defined in 
exactly the same way within the basis of domain wall 
functions (defined in Eq. iJTJJ), with the domain wall 
wavevectors k an (Eq. I|18a|) ') in place of k an . 

The total scattering matrices for coherent propagation 
through the disordered region in the case of a domain wall 
and uniform magnetization can be written respectively as 



4o m) = Pl*Si*p 2 *---*S Ni * p Ni + l- 



(12a) 
(12b) 



For the incoherent conductance, the phase coherence 
between successive scattering events is assumed to be 
lost. The resulting scattering matrices for the domain 
wall and uniform cases are obtained from the combina- 
tion law of Eq. ©, but using probabilities (H) instead of 
amplitudes (£). We then write 



-r(dw) 



S (L) *S!*S 2 *---*S Ni * S (R) , (13a) 
*S Ni . (13b) 



s£f = s 1 *s 2 * 
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The matrices p a and p a do not appear in this case since 
they contain a pure phase shift and the associated trans- 
mission probabilities are simply unity. 

In the quasi-ballistic regime, the elastic mean free 
path is much larger than the system size and only first- 
order scattering processes are important. For the uni- 
formly magnetized case the incoherent and coherent con- 
ductances become equal in this regime. In the diffu- 
sive regime, on the other hand, constructive interference 
of time-reversed paths with identical starting and end 
points leads to an enhancement of the coherent reflection 
known as weak localization^ In the limit of large disor- 
der, the coherent conductance enters the strongly local- 
ized regime, where g decreases exponentially with length 
and the system becomes insulating^ By constrast, the 
incoherent conductance scales like the classical Drude 
conductance (g ~ 1/L), even for arbitrarily large dis- 
order. 

A characteristic feature of phase coherence effects in 
the diffusive regime (such as weak localization or conduc- 
tance fluctuations) is their universality, or independence 
of system size, which means that their relative impor- 
tance is largest for small system sizes. In our model this 
is particularly significant since we are constrained to work 
with system sizes that are smaller than typical systems on 
which experiments are performed. Therefore, the relative 
importance of coherence effects appears overemphasized. 
On the other hand, the incoherent conductance scales 
with the system size and our numerical results can be 
safely extrapolated to larger sizes. Furthermore, in the 
incoherent case the fluctuations of conductance decrease 
with increasing disorder, which allows average quantities 
to be computed accurately with relatively small numbers 
of samples. 

A fundamental limitation with finite numerical mod- 
els of coherent transport is that, for a given width L y 
and impurity density, there is a maximum value of the 
system length L z beyond which the system enters the 
strongly localized regime. However, if L z is larger than 
the phase coherence length L$ then the localized be- 
haviour obtained in the coherent model is not relevant. A 
comparison between coherent and incoherent results per- 
mits us to determine those features which are due to co- 
herence and will disappear with increasing temperature. 
At non-zero temperature the presence of phase-breaking 
scattering means that the relevant regime for transport 
is intermediate between the coherent and incoherent lim- 
its. It is in general very difficult to treat this regime, 
but by considering both coherent and incoherent limits 
we are able to gain some insight into the behaviour of 
experimentally relevant cases. 



B. Electronic states in a ballistic domain wall 

For our model of a domain wall we assume a linear 
rotation of f(r) = f(z) in the yz plane over the region 



< z < A, described by 

f (0, sin[0(»], -cos[0O)l), 0<z<A, 
f{z) = \ ~ ~ (14) 

I (0, 0, sgn[z]), otherwise, 

where 6{z) = ttz/X is the angle between the magnetiza- 
tion inside the wall and the z-axis. Using this profile, one 
can describe the qualitative features of transport through 
a domain wall, which have been shown to be independent 
of the detailed form of the wall. 15 In addition, the profile 
of Eq. H14f) has the advantage that the corresponding ba- 
sis states can be found in closed form. 15,1 — The local spin 
eigenstates |±) (parallel and anti-parallel to the direction 
of f(z)) are position-dependent inside the domain wall. 
Denoting the spin basis states with respect to the fixed 
z-axis as |±) , we have 



where 



R[6(z)] = 



= R[6(z)] 



cos[0(z)/2] 
sin[0(z)/2] 



-sin[0(z)/2]\ 
cos[0(z)/2W 



(15) 



(16) 



is the spinor rotation operator containing the amplitudes 
{u'\a) representing the transformation from fixed to lo- 
cal basis. Because the states |±) are position-dependent 
inside the wall, the eigenfunctions of Eq. Q in this re- 
gion are not spin eigenstates, but a combination of both 
up and down components: 



v4(r) = 



mv 



s ct> n (y)(\<T)±iA an \-a)). (17) 



Here the symbols ^ denote the direction of electron prop- 
agation (i.e. right or left), while a = ± is a quantum 
number representing the spin of the state ip a n (r) to which 
ipa-n('r) reduces in the limit A — > oo. The other parame- 
ters are 



ha 
A a 



k 2 



16fc^ 



2fc>fc, 



k 2 



k 2 . 



hk„ 



i + Ai 



(18a) 
(18b) 

(18c) 



with k 0n = y (fc+.„ + k 2 _ n )/2 and k x = n/2X. Note that 
the relation between the velocity factor v an and wavevec- 
tor kan is more complicated than for the uniform basis 
states in Eq. (JSJ). 

The number of propagating modes per spin quantum 
number inside the domain wall region, N a , is given by 
the maximum value of n for which k an is real- valued. For 
parameter values of interest we generally have N a = N a . 
However, for very small A the rotating potential can lead 
to effective band gaps (analogous to those occurring in 
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spin-independent periodic potentials), leading to N a < 
N a . An additional complication for small A is that states 
with complex k an , which are inconvenient to treat within 
our scattering matrix approach, may become relevant 4£ 
Fortunately, such small values of A do not fall within 
the range of parameter values for the systems we are 
considering, and so in our calculations we always have 
N a = N a with no complex wavevectors. 

In the domain wall model defined by Eq. 114|) . elec- 
trons incident from the left at z = or from the right at 
z = A are scattered due to the change from uniform to ro- 
tating magnetization. The scattering matrices for these 
interfaces, which we write and , can be calcu- 
lated using the standard method of matching incident 
and scattered wavefunction components. Since the do- 
main wall potential is uniform in the transverse direction 
(i.e. has no y dependence), the interface does not mix 
different transverse modes. We thus have = if 

n' 7^ n, for £ = t, r, t' ', r'. The diagonal amplitudes £^ R ^„ 
are determined by calculating the appropriate scattering 
state solutions for states incident on the interface, lead- 
ing to 4 x 4 sets of linear equations which we present in 
Appendix El (Eqs. QAlHA2[) j. 

For general parameter values it is most convenient to 
solve Eqs. I|A1IIA2|) numerically. However, simple asymp- 
totic expansions can be found in the limits of wide 3 ^ and 
narrow — walls which provide useful insight. In the exper- 
imentally relevant case of small spin-splitting and large A, 
a particularly simple solution can be found by expand- 
ing to second order in the inverse of the dimensionless 
parameters 



Pn = k\/2k 0n kx, 
Pf = k%/2k F k\. 



(19a) 
(19b) 



These parameters characterize the effective "width" of 
the domain wall for an electron in transverse channel n 
(p n ) and for the wire as a whole (pf)- For an individual 
channel n, the adiabatic limit corresponds to p n 3> 1 and 
is most readily obtained for channels with large n (and 
hence small fcon)- As discussed in Ref. [TEL the degree 
of adiabaticity is in general channel-dependent. Never- 
theless, pf permits a characterization of adiabaticity for 
all states: for a wide wall we have pp > 1, and hence 
p n > 1 for all channels. Within this assumption we solve 
Eqs. (|A1IIA2[) to 0(l/pD for each n, obtaining 



~(L,R) = ,(L,R) 
cr'n' ;an a'n' ' \crn ' 



(20a) 



~(L,R) ~/(L R) g L_ J_\ (2Qb) 

L cm';an ~ l an':an °« " I 1 J ' ^ UL 7 

f (L3) - ~t' (L ' R) -=S , — (2Qc) 

b -an':crn ~ b -an';<Tn h°n'n ■ ^ZUCJ 

Pn 

Eqs. H20|) show that there is essentially no reflection for 
electrons incident on the domain wall, but that the main 
effect of the interfaces is to scatter electrons into a su- 
perposition of up and down transmitted channels (con- 



serving the mode number), with amplitude determined 
by p n - 

In the regime of small splitting, A <C Ep, the trans- 
port properties of the domain wall are determined pri- 
marily by the parameter pp, while the dependence on 
E-p, A, and A individually can, to a good approximation, 
be neglected. This is true both for the intrinsic domain 
wall scattering shown in Eqs. I|2U|) and for the impurity 
scattering to be discussed below. For large splitting, 
the dependence on A becomes important as the differ- 
ence in the number of up and down conducting channels, 
N + — becomes significant. In this work, however, we 
are interested only in the case of weak splitting which is 
relevant for transition metal ferromagnets. 



1. Coherent and incoherent conductance for a ballistic wall 

Before introducing disorder in the following section, we 
now briefly consider the case of a disorder-free ballistic 
domain wall (V(r) = 0), which can be treated by setting 
Ni = in Eq. (|12a|l (coherent case) or Eq. H13a(l (inco- 
herent case). An incoherent ballistic system represents 
an idealized scenario in which electrons are subject to 
phase-breaking events which do not affect momentum. 33 
As such, it is a useful way to treat the effects of deco- 
herence in mesoscopic transport. In the case of a do- 
main wall, it is important to understand the differences 
between incoherent and coherent transmission in the bal- 
listic this gives insight into effects which are also 
relevant for the disordered case. 

Within the asymptotic approximation of Eqs. J3UJ, the 
transmission probabilities for coherent and incoherent 
ballistic walls are, to 0(l/p n ), given by 




Sn'nS- sin 2 , coherent, 

(21a) 



incoherent, 



1 an' :tm 



crn 1 :crn • 



(21b) 



From Eqs. 121|) the basic difference between the coherent 
and incoherent cases is a suppression of the oscillatory 
component of the spin-dependent transmission in the lat- 
ter case. This occurs because the oscillations observed in 
the coherent case arise from the phase interference be- 
tween up and down wall basis state components com- 
prising the electron scattering state, which is suppressed 
in the incoherent case. This result implies that the oscil- 
latory torques exerted by a s pin- polarized current on a 
domain wall predicted in Ref. |32J would be suppressed if 
the transport were incoherent. 

In Figure El we show the exact total spin-dependent 

transmission probabilities T ++ and T |_ (defined in 

Eq. (|8a|l ) as a function of p?, for both coherent and in- 
coherent ballistic walls. The agreement with the asymp- 
totic results from Eqs. 112 1H is shown in the inset; it can 
be seen that these are accurate for moderate to large pp 
values, but for small p-p they diverge significantly. In 
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FIG. 2: Spin-dependent transmission per incoming channel, 
T a i a /N a , for a ballistic domain wall as a function of the di- 
mensionless effective wall width pp (bottom axis). For com- 
parison, the corresponding actual wall width A for parameter 
values Ef = lOeV and A = O.leV is also shown (top axis). 
The curves presented are for a — + , with a' — ±<r, and are 
indistinguishable from the corresponding ones for a — — (for 
the parameter values used). Solid and dashed lines indicate, 
respectively, coherent and incoherent combination of scatter- 
ing matrices. The actual wire width used is L y = 5nm, cor- 
responding to N± = 25, but the curves are indistinguishable 
from their form in the limit of large L y . Inset: comparison of 
the asymptotic solutions for T aa from Eq. (I21bt (thick lines) 
with the exact values (thin lines) over a larger range of j>f- 



the limit pf — ► 0, the exact solutions show that the 
incoherent transmission goes to N a /2, which is some- 
what surprising since we expect complete mistracking 
(i.e. T aa = 0) in this limit, as is observed for the coherent 
case. The origin of this result can be understood as fol- 
lows. By combining scattering matrices incoherently we 
are effectively "measuring" the individual path taken by 
an electron going through the wall, so that is projected 
onto one of the right-moving wall basis states ip±(^)- 
Since each of these states comprises an equal weighting 
of local up and down components as pv — ► 0, the incoher- 
ent combination must also yield an equal weighting. It 
is clear that the incoherent result is unphysical for small 
Pf, since there will always be some non-zero distance 
over which transport is coherent. We therefore need to 
be careful when interpreting the incoherent results for 
small pf ■ 



C. Delta function model of scattering 

In transition metal ferromagnets the elastic spin- 
conserving scattering of conduction electrons arises from 
a variety of mechanisms, including impurities, defects, 
and grain boundaries. These scattering rates depend 
on the underlying details of band structure, and are in 
general strongly spin-dependent. In spite of this com- 



plex situation, the strength of elastic scattering can be 
characterized by only two parameters, namely, the spin- 
dependent elastic mean free paths l±. The precise mi- 
croscopic model leading to a given value of l± is then 
of minor importance. In the sd model the spin depen- 
dence of l a arises from two sources^ the intrinsic spin 
dependence of the electron wavevectors k an and a differ- 
ence in the number of available d states of each spin into 
which the conducting s electrons can be scattered. The 
former effect is due to the spin-splitting of the s band, 
which makes fc +i „ > fc_ )n . Therefore, electrons with spin 
up are generally less strongly scattered than those with 
spin down, which in the absence of other sources of spin 
dependence leads to Z+ > Z_. However, since we work 
in the regime of small splitting, A <C Ep, this difference 
is rather small. The dominant contribution then comes 
from scattering into the d band. The resulting spin de- 
pendence can lead to Z+ > Z_ or Z + < Z_ depending on 
the form of the up and down d sub-bands at the Fermi 
energy. For definiteness, we assume that the spin down 
d sub-band has a greater number of states at Ef than 
the spin up sub-band, which implies that spin down elec- 
trons are the more strongly scattered ones, leading to 
l + > i_. To represent scattering processes, we adopt a 
simplified picture in which s electrons are scattered by 
the static potential V(r). The degrees of freedom cor- 
responding to the d electrons are therefore not explicitly 
included and, in particular, we ignore the possibility that 
electrons scattered into d states might not return to the 
s band. 

The potential V(r) is defined as a random array 
of delta functions with spin-dependent amplitudes u±, 
where the orientations are defined with respect the lo- 
cal magnetization direction f(z). The delta function 
model is a convenient phenomenological approach to im- 
purity scattering which has been widely used in spin- 
independent mesoscopic transport theories^iS as well 
as in ferromagnetic systems i&SH The spin dependence 
of scattering in the ferromagnet is determined by the 
ratio of up and down amplitudes, which we write as 
p = u-/u + . For spin-dependent impurity scattering we 
take p > 1, which corresponds to > Z_. For complete- 
ness, we will also consider the case of spin-independent 
impurity scattering, p = 1. 

We write the up and down components of V(r) in the 
local spin basis as 

Ni 

V a (r) = Ucr ^2S(?-r a ), (22) 

a=l 

where there are N- x impurities with positions r a randomly 
distributed in a region of area L y xL z . In general, L z may 
be different from the domain wall length A. The impurity 
density is n\ — N{/L y L z . The total impurity potential 
can then be written as a sum of spin-independent and 
dependent terms: 

V+(r) + V (r) Vj-(r) - V (fO - 
V{r) = V+[ )+ 2 [ > t + V+[r> 2 [ ) f(z) ■ S. (23) 
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Alternatively, the potential can be written in a diagonal 
form 

V(r) = iT 1 [*(*)] ■ v ° {7 ^ ■ R[0(z% (24) 

in terms of the rotation matrix R[Q(z)\ defined in 
Eq. P|. 

As shown in Eqs. l|12H13f) . the transmission through 
the disorder potential can be calculated by combining 
the scattering matrices s a , s a of the individual delta 
function scatterers located at f* Q , together with the ma- 
trices p a and p a for propagation between successive scat- 
terers. The dependence on the longitudinal component 
of the scatterer position, z a , is contained in p a and p a , 
while the dependence on the transverse component, y a , 
is contained in s a and s Q . We therefore write the scat- 
tering matrix for a delta function scatterer located at f a 
as s^(y a ) (uniform) and s^(y a ) (domain wall). We 
describe the calculation of these scattering matrices in 
Appendix 151 and present solutions valid in the Born ap- 
proximation in Appendix [HI 

To understand the scattering from a delta function 
with spin-dependent amplitudes u a , it is useful to con- 
sider the total spin-dependent reflection for the uniform 
case, R&}. Figure EH shows (R&}), the reflection averaged 
over transverse positions y a of the scatterer, as a func- 
tion of u,y for several values of p. For small u ai we see 
that (R [ J}} goes as u 2 ,, which can also be seen directly 
from the Born approximation^ The ratio of up to down 
back-scattering thus goes as p 2 for small u a . For larger 
u a , however, (R&}) is bounded above by 1, attaining this 
value in the limit u„ — ► oo. This occurs because a delta 
function scatterer can "block" at most one channel^ 
The spin dependence of reflection is then dramatically 
reduced for large u a . Since our goal in this work is a 
model with spin-dependent scattering, we must therefore 
work in the regime of small u a . However, we cannot take 
u a arbitrarily small, since decreasing u a means that a 
larger N a is required to achieve a given disorder density, 
leading to an increase in computation time. In practice 
we have found u + < to be a useful compromise for 

parameters in the range 1 < p < 2. 



III. TRANSPORT THROUGH A DISORDERED 
DOMAIN WALL: INTRINSIC PROPERTIES 

In this section we consider the intrinsic conductance 
properties of a disordered domain wall. Here the domain 
wall constitutes the entire system of interest, i.e. L z = A, 
and is therefore assumed to be connected on either side to 
perfect (disorder-free) leads. By comparing with a uni- 
formly magnetized disordered region of the same length, 
we calculate an intrinsic domain wall magnetoconduc- 
tance which measures the difference in conductance due 
to the presence of a domain wall. We consider both coher- 
ent and incoherent regimes; the corresponding scattering 
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FIG. 3: Total spin-dependent reflection from a delta function 
scatterer with spin-dependent amplitudes u a in the uniformly 
magnetized case, (_rS) (averaged over transverse positions 
y a ). Curves are shown as a function of (u+/Ef) 2 for different 
values of the up/down scattering ratio p = u~/u+. For p = 1 
the only spin dependence of scattering is due to the small 
difference between k+, n and k- >n , and hence — {R+\_}. 

For small u a , goes as v% (Eq. IC4al h while for large u a 

it approaches a maximum value of 1. The particular system 
used in this calculation has L y = 2nm, giving iV± = 10. 

matrices for the system are obtained from Eqs. i|12H13[l . 
Our primary aim in this section is to study the role of 
impurity scattering inside the domain wall. The case of 
most interest is spin-dependent disorder (p > 1), how- 
ever, in order to relate our model to previous works and 
more clearly understand aspects of the spin-dependent 
case, we also consider spin-independent disorder, p = 1. 

In the uniform case, the two spin sub-bands of con- 
duction electrons are uncoupled, so there is no spin 
mixing. All off-diagonal (in spin quantum numbers) 
transmission and reflection amplitudes are therefore zero, 
i.e. £- CTn ';o-n — for £ = t,r,t',r'. The total spin- 
dependent transmission and reflection probabilities (de- 
fined in Eqs. JHJl ) satisfy 

S (uni) =0j S (uni) =s (uni) 5 (2ga) 

for Z = T,R, T', R'. 

The conductance properties of the up and down sub- 
bands of the uniform system are well described by ran- 
dom matrix theoryj2£ii2& whose expressions for the average 
transmission and its moments are in excellent quantita- 
tive agreement with the calculations performed with our 
model. For a system of length L z , the total transmission 
for spin sub-band a in the incoherent case is given by 

(Ti2>=JV(l + Vk) s (26) 

where N a is the number of channels and l a is the spin- 
dependent elastic mean free path. Eq. (|2*rj|l expresses 
the transmission as a series combination of the ballis- 
tic contact transmission, N„, with the Drude conduc- 
tance, N a l a jL z . In the diffusive regime, l a <C L Zl the 
Drude term dominates and (T^™^) displays an Ohmic 
1/L Z dependence. In the quasi-ballistic regime, l a 3> L z , 
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Eq. (|26|l is only approximately correct, but is within 
several percent of the value obtained by a more precise 
calculation^ 

For delta function scatterers in a two-dimensional 
quasi- ID geometry, the mean free path l a appearing in 
Eq. l(2T)| is given by 



where w CT! f = Kh^y/m is the Fermi velocity in spin sub- 
band a. This definition is a factor of 2 larger than the 
two-dimensional formr 8 since in a quasi-lD geometry for- 
ward scattering processes do not reduce transmission^ 

The mean free path l a can be interpreted intuitively 
as the typical distance travelled by an electron of spin a 
before undergoing a momentum randomizing scattering 
event. However, it should be noted that such a concep- 
tual scattering event does not correspond to the scatter- 
ing from an individual delta function scatterer in our 
model: the length l„ typically corresponds to a large 
number of individual scatterers. 

The coherent transmission, (T±^j ), is reduced with 

respect to (T±fa C ) by the weak localization correction, 
which in a quasi-lD system has the limiting value —1/3 
in the diffusive regime L z /l a ^> 

A. Disorder-induced enhancement of 
spin-mistracking in domain wall transport 
coefficients 

We first consider how the transport coefficients of a 
domain wall vary as a function of disorder. It is most 
instructive to study the total spin-dependent transmis- 
sion and reflection, summed over all transverse channels 
and averaged over impurity configurations, (T^f™ ) and 

(i?^). In Figure 01 these quantities are shown for a rel- 
atively wide wall (p-p = 5) with spin-dependent disorder 
(p = 2) as a function of disorder strength (1/1+ in units 
of 1/A). Both coherent (Eq. (|12a[) ) and incoherent cases 
(Eq. Ijl^al) - ) are shown. 

At zero disorder, 1/1+ = 0, transport through the bal- 
listic wall is highly adiabatic, so that (Taf^) 3> (T^™J). 
With increasing disorder, the dominant transmission co- 
efficient, (tJct W ^), decreases quite rapidly, with a form 
similar to that of the transmission in the uniform case, 
(Tjr } ) (Eq. HI). Also, we observe that (Tl d _ w) ) de- 
creases more rapidly than (T++ }, since Z_ < l+. These 
initial decays can be quantitatively described by the Born 
approximation of Appendix [U] and the coherent and in- 
coherent cases yield very similar results. 

The off-diagonal (in spin) coefficients (T^J) also ex- 
hibit a linear decrease with 1/1+ in the quasi-ballistic 
regime. However, for a relatively weak disorder (X/1+ — 
1) the behaviour of off-diagonal coefficients begins to dif- 
fer considerably from that of the diagonal ones. As we 



can see in the inset of the upper graph in Figure 0] the 
negative slope of the coherent (T_J$) levels off, and the 

incoherent (T_J*J) increases with disorder. 

In both cases (coherent and incoherent), the magni- 
tude of (T^'J) remains relatively constant as a function 
of disorder. This means that the relative transmission 
with spin-mistracking, (T^J)/({T^) + {T<%})), in- 
creases dramatically as a function of disorder. Further- 
more, since (T^) < (T} d + w) ) while (T^) = (Tl d + w) ), 
the proportion of transmission with mistracking is greater 
for the spin down than for the spin up sub-band. 

The case of spin-independent disorder, p = 1, (not 
shown) is qualitatively similar to the spin-dependent case 
just discussed, with the difference that (T++) ~ (T__) 
and (R++) — (R--). This is because in this case the 
only spin dependence of the scattering arises from the 
small spin dependence of the wavevectors k an . 

The preceding observations imply that the adiabatic- 
ity which applies to ballistic domain wall transport no 
longer applies in the presence of disorder. In particular, 
a wall which is highly adiabatic in the regime of ballistic 
transport becomes less so in the presence of disorder. For 
spin-independent disorder, a similar result was recently 
obtained in the context of transport through disordered 
wires in the presence of inhomogeneous magnetic fields. 40 
For spin-dependent disorder, the reduction in adiabatic- 
ity depends on the spin direction, so the wall can no 
longer be characterized by a single adiabaticity param- 
eter as it was in the ballistic case^ The origin of the 
enhancement of mistracking with increasing disorder can 
be understood intuitively as the cumulative effect of small 
amounts of mistracking acquired in scattering from each 
individual delta function (see Eqs. I | (J2b .dll. We illus- 
trate this important idea with a simple one-dimensional 
toy model in Appendix ID1 

Finally, we note several differences between the coher- 
ent and incoherent cases in Figure^ In the ballistic case, 
1/1+ = 0, there is a small difference between (T^^ co ) and 

{Tpf™l nc ). This is due to suppression of the oscillatory 
component of the coherent transmission, which was dis- 
cussed in Section III B II and illustrated in Figure [21 This 
difference may be positive or negative, depending on pp, 
and becomes increasingly significant as pp — > 0. The 
effect persists for small but disappears for larger 

disorder since the phase information corresponding to 
the precessional component is lost after many scattering 
events. 

A second difference between the coherent and incoher- 
ent cases is that for large disorder the scattering coeffi- 
cients with mistracking, (T^fJ) and (R^™}), are smaller 
for the coherent case than for the incoherent one. The 
precise origin of this difference is not clear, but we can 
eliminate several possible reasons. Firstly, by looking at 
the equivalent curves for different values of p-p , it is found 
that the difference has a constant sign for all A. This sug- 
gests that it cannot be explained by the suppression of 
the precessional component of transmission in the inco- 
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FIG. 4: Spin-dependent domain wall transmission and reflec- 
tion, (Tain) and (ii CT / CT ), (normalized by the number of con- 
ducting modes per incident spin channel, N a ) as a function of 
disorder, measured by 1/1+ (shown in units of 1/A). Coher- 
ent (solid lines) and incoherent (dashed lines) cases are shown 
(the arrows connect corresponding values) . The wall width is 
Pf — 5, which corresponds to a wide wall close to the adiabatic 
limit {T aa > T-aa for the ballistic wall at 1/1+ = 0). The pa- 
rameters for the delta function scatterers are u+ = Ef/5 and 
p — 2, and the averages are performed using Ns = 2000 im- 
purity configurations. The wire width is L y = lOnm, leading 
to N± = 51. 



herent invoked in the discussion of the previous 

paragraph. Furthermore, it is found that the difference 
scales linearly with the number of channels N a (or equiv- 
alently L y ). It cannot, therefore, be explained as a weak 
localization effect, which should be characterized by a 
constant magnitude, independent of L y . 



B. Intrinsic domain wall magnetoconductance: 
sign reversal in coherent case 

We define the intrinsic domain wall magnetoconduc- 
tance as the difference between the conductance of a uni- 
formly magnetized region ( "uni" ) and that of a domain 
wall ( "dw" ) of equal length with identical impurity con- 
figurations 



A.g = g u 



Sdv 



(28) 
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Assuming that we can contact the wall directly, the im- 
purity average (Ag) could be associated with the change 



FIG. 5: Coherent (solid lines) and incoherent (dashed lines) 
domain wall magnetoconductance (Ag) (filled symbols) as a 
function of 1/1+ (in units of 1/A) for a domain wall with 
Pf = 5, in the case of (a) spin-independent disorder (p = 1) 
and (b) spin-dependent disorder (p = 2). The differences 
in spin-dependent transmission (T+) are also shown (empty 
symbols) . The system parameters are as for Figure 2] with 
iVs = 3000 impurity configurations. 



in the measured conductance when an external magnetic 
field along the z direction is applied in such a way as 
to destroy the domain wall and arrive at a magnetically 
homogeneous configuration. Although they are not of di- 
rect experimental relevance, we also define the differences 
in spin-dependent transmission, AT CT = j"j um ) — Ti dw ^ in 
order to guide our physical discussion. With the above 
notations we obviously have Ag = AT+ + AT_ . 

In Figure[S]we show the disorder-averaged (Ag) (thick 
lines, filled symbols) and (AT±) (thin lines, empty sym- 
bols) as a function of disorder, for spin-independent 
(p = 1) and spin-dependent (p = 2) disorder, in both 
coherent and incoherent cases. 

In the spin-independent case (p = 1, Figure [S^,), the 
main feature is a negative coherent magnetoconductance, 
which becomes positive and very small in the incoher- 
ent regime. For the spin-dependent case (p = 2, Figure 
\^jp) the negative coherent magnetoconductance is obtained 
above a threshold disorder, and a positive (Ag co ) appears 
in the incoherent regime. Below we comment on the gen- 
erality of these basic findings, their relationship with pre- 
viously found effects, and their physical relevance. 

A negative coherent magnetoconductance was pre- 
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dieted by Tatara and Fukuyama-2«ii for the case of spin- 
independent disorder, and interpreted as a weak localiza- 
tion effect. Such an effect has, however, eluded experi- 
mental confirmation. The underlying reason for the pu- 
tative reduction of weak localization in a domain wall is 
a suppression of backscattering processes which conserve 
spin. That is, r an i. an is reduced due to the possibility of 
scattering into the opposite spin channel, r_ CTn ' ;o - n . Since 
weak localization stems from an enhancement of the di- 
agonal spin-conserving reflection amplitudes r an . an (co- 
herent backscattering), it follows that the effect will be 
reduced in a domain wall as compared to a uniformly 
magnetized region. 

Our numerical results suggest that a suppression of 
weak localization by the domain wall is indeed the dom- 
inant mechanism responsible for the coherent magneto- 
conductance in the regime of large disorder. Firstly, the 
limiting value (Ag co ) ~ —0.1, obtained in the diffusive 
regime, is of the order of the quasi-lD weak localization 
value of —1/3 for the uniform caseiSS^ Furthermore, this 
limiting value is approximately independent of system 
size (L y and \)r^ which is a general characteristic of 
coherence effects such as weak localization. 

From the point of view of physically measurable effects, 
it is important to note that the coherent magnetoconduc- 
tance (Ag co ) is characterized by relatively large fluctu- 
ations. Indeed, in the diffusive regime Ag co follows ap- 
proximately a normal distribution with fluctuations char- 
acterized by ^Var(A 3co ) = ((A 5c 2 Q ) - <A 5co ) 2 ) 1 /2 ~ .3, 
as illustrated in Figure EJ For large disorder, the magni- 
tude of these fluctuations is independent of system size 
{L y ) and up/down scattering ratio (p), and is related to 
the universal fluctuations of <?dw and g U ni- We discuss 
these aspects in more det ail in the following section. We 
notice that ^/V&r(Ag co ) ~ 3|(A<? co )| for large disorder. 
For an individual disorder configuration, there is thus 
a significant probability for A<7 co to be positive as well 
as negative. This could make it difficult for a negative 
(Ag co ) to be detected experimentally (even in the coher- 
ent regime) beyond statistical uncertainty. 

Another significant factor concerning the coherent 
magnetoconductance is that in the case of spin- 
dependent disorder (Ag co ) is positive for small disorder. 
The value of 1/1+ at which it changes sign increases with 
the system width L y . In our calculations, numerical con- 
straints limit us to systems containing on the order of 10 2 
channels. However, the nanowires of experiments such as 
Ref. contain on the order of 10 4 channels. The region 
in which (Ag co ) is negative would then correspond to an 
unrealistically large disorder, particularly for the cobalt 
nanowires of Ref. 4, in which domain walls are relatively 
narrow (A ~ 15nm). For materials with wider walls, such 
as Nickel (A ~ lOOnm), the necessary density of impurity 
scatterers would be smaller, although it would be more 
difficult to have phase coherence across the greater length 
of the wall. Thus, to access the regime of negative (Ag co ) 
experimentally it would be necessary to work with very 
small, highly disordered nanowires at low temperature. 
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FIG. 6: Distribution of the coherent intrinsic magnetocon- 
ductance, Ag co , for spin-independent (p = 1, x) and spin- 
dependent (p — 2, +) disorder. The curves in the cen- 
tre show (Agco), while the upper and lower curves show 
(Agco) ± Va,r(Ag co ) 1 ^ 2 respectively. All parameters are as 
for Figure 



In the incoherent regime, the magnetoconductance 
(Ag- mc ) in the case of spin-independent disorder is pos- 
itive but extremely small (approximately two orders of 
magnitude smaller than KAT-^mc)!). For spin-dependent 
disorder (A<?; nc ) is much larger. Although not shown in 
Figure El we note that the incoherent quantities scale 
linearly with increasing L y , while for fixed 1/1+ they de- 
crease with increasing A (see Figure EJ- 

The behaviour of (Ag- mc ) can be understood as a com- 
bination of spin-mistracking with spin-dependent scat- 
tering. As we saw in Section IIII Al successive scatter- 
ing events in the domain wall enhance spin-mistracking. 
This means that the electron is sensitive not only to 
the impurity potential associated with its incoming spin, 
but also to the one with the opposite orientation. Since 
the scattering strengths are spin-dependent, this leads 
either to a reduction or enhancement of transmission 
compared to the uniform case, according to the spin di- 
rection: (A7+, inc ) > and (AT_, inc ) < 0. For spin- 
independent disorder, these differences are small and al- 
most equal in magnitude, so that (Agi nc ) is negligibly 
small. For spin-dependent disorder, however, we have 
|(AT +jinc )| > |(AT_ iinc )|, which leads to a positive and 
much larger value of (A<?j nc ). In this case there is thus a 
significant magnetoconductance arising from an enhance- 
ment of back- scattering due to the exposure of electrons 
in a superposition of spin states to spin- dependent scat- 
tering. This idea was used at a phenomenological level in 
Ref. IliL where transport through a diffusive domain wall 
in the almost adiabatic limit of small spin-mistracking 
was described by a reduced effective mean free path rep- 
resenting a weighted average between spin up and down 
mean free paths. 

In the ballistic limit, 1/1+ — > 0, there is no impurity 
scattering inside the wall. The back-scattering is then 
entirely due to reflection from the wall interfaces, which 
is extremely weak. Both (Ag co ) and (Agi nc ) therefore 
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become negligible in this limit. 

In the diffusive regime (X/1+ <C 1), we believe that two 
primary mechanisms explain the different behaviour of 
(Ag co ) and (A<?i nc ). Firstly, the limiting value of (Ag co ) 
for large disorder appears to be due to a reduction of 
weak localization by the domain wall, which does not 
apply in the incoherent case. Secondly, as we discussed 
in Section IIII Al the spin-mistracking transmission and 
reflection, (T^J co ) and (R { ^ co ), arc reduced in the 
coherent regime with respect to the incoherent one. We 
propose that this leads to a suppression of the enhanced 
backscattcring in the domain wall which is dominant in 
the incoherent case. For large disorder, the positive mag- 
netoconductance present in the incoherent case is then 
cancelled out in the coherent case, leaving only the neg- 
ative component from the weak localization reduction. 

C. Reduction of universal conductance fluctuations 
for coherent domain wall 

An important feature of quantum transport through 
diffusive coherent systems is that the conductance fluctu- 
ations, (Sg 2 ) 1 / 2 = Var(g) 1 / 2 , are universal, with a mag- 
nitude of order e 2 /h independent of the system size or 
mean free path4^ For a spinless quasi- ID disordered sys- 
tem, the conductance fluctuations have the value y/2/15 
in units of e 2 /h. In our model this applies to each spin 
sub-band in the uniformly magnetized case, so that 

(5[T£ uni) ] 2 ) 1/2 = \/27T5- (29) 

For a non-magnetic system ("rim"), there is no spin- 
splitting of the s band (which would correspond to A = 
in our model) and the transport of up and down elec- 
trons is identical. In that case we have Tjl™^ = yt 11111 ) 

and hence {6[T<? m) + T^} 2 ) 1 / 2 = 2{5[T + nm) ] 2 ) 1 ' 2 . The 
conductance fluctuations are then 

(^ m ) 1/2 = mT [ r ] \ 2 ) 1/2 - (so) 

In a ferromagnetic system (A ^ 0), on the other hand, 
the spin dependence of the wavevectors k an suppresses 
the correlation between j'j" 111 ) anc l j^ unl ) _ j n ^his case we 



have (a[Tj^+r< uni) ] 2 )V2 = {25[T + m{) ] 2 f/ 2 , so the con- 
ductance fluctuations are reduced as compared to those 
of a non-magnetic system: 

(^L) 1/2 = (2^ ) ] 2 ) 1 / 2 = V47l5. (31) 

This reduction due to the non-degeneracy of spin states is 
directly analogous to the experimentally observed reduc- 
tion of conductance fluctuations due to Zeeman splitting 
in an applied magnetic field^ 

In FigureJTJi we show the conductance fluctuations ob- 
tained from our model in the uniform case, {Sg 2 ^^) 1 ' 2 and 

(S[T+ uni) ] 2 ) 1 / 2 , as a function of 1/1+. They arc in excel- 
lent agreement with the theoretical values from Eqs. I|29|) 
and for X/1+ > 2. 

Figure also shows the conductance fluctuations 
for a domain wall region, (Sg^) 1 ^ 2 , together with 
(j[ T | dw )]2)i/a. We see that the conductance fluctuations 
in a disordered domain wall are reduced with respect to 
those of a uniformly magnetized region. Moreover, the 
conductance fluctuations are no longer universal since a 
slow decrease is observed as a function of 1/1+ in the 
diffusive regime. 

This result can be understood as arising from statisti- 
cal decorrelation between the components of J^ dw ) . For 
the domain wall coefficients we do not have a way of 
estimating the fluctuations analogous to Eq. J2HJ. How- 
ever, we can make the hypothesis that the relative fluctu- 
ations for each component are the same as in the uniform 
case: 

(5[T^} 2 ) 1/2 „ {5[T^] 2 f/ 2 

(T (dw) ) (T (uni) ) 

Furthermore, we can expect that T^"^ and T^™J will be 
uncorrelated in the diffusive regime, so that 

( tf[T< (dw)]2 )1 /2 „ (S[T£^} 2 + 8[T { ^] 2 ) 1 ' 2 . (33) 

Using the fact that (T^ dw) ) = (t£ w) ) + (T^J), Eqs. G2h 
I33|l then give 



((5[T( (dw) ] 2 ) l/2 



(dw), 

I / 

(uni). 



9/ 7 i( d »)\/ T ( dw )l\ 1 ^ 2 
~ ' " ' -n-n-/ \ ^jr(uni)j2\l/2 



(dw), 2 



(34) 



The first factor on the r.h.s. is approximately unity since 
(Ti dw) ) ~ (Ti uni) ). We thus see that (S[T^ dw) ] 2 y/ 2 is 
reduced with respect to (JfTi" 111 '] 2 ) 1 / 2 by the second term 
on the r.h.s. of Eq. (|34|l . This reduction factor is most 



important when the two transmission coefficients (T^o- ) 
and (T^fJ) are of the same order, which is the situation 
approached with increasing disorder. 

In Figure [7|d we compare the actual fluctuations 
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FIG. 7: (a) Fluctuations of conductance, (Sg 2 ) 1 / 2 , (solid 
lines) and spin-dependent transmission, (5T+) 1 / 2 , (dashed 
lines) for a uniformly magnetized region (o) and domain 
wall (A), in the case of spin-independent disorder (p — 1). 
The fluctuations of spin-down transmission, (ST 2 ) 1 / 2 are not 
shown, but agree with (5T+) 1 / 2 to within the statistical error 
of the calculation. The dotted horizontal lines indicate the 
theoretical values in the uniform case, (Sg 2 ni ) ll/2 = y/ 4/15 

and (<5[T± 11 '] 2 ) 1/ ' 2 = \/2/15, which are in good agreement 
with our numerical results for X/1+ > 2. (b) Fluctuations 
of the spin-dependent domain wall transmission coefficients: 

(5[r |dw) ]2)1/2 (A)) (5[r |dw) ]2)1/2 ( + ^ Wr (dw) ]2)1/2 (>< ^ 

The corresponding estimated values, based on the approxi- 
mation in Eqs. (I32H34H . are also shown (solid lines). System 
parameters are pf — 5 and L v = lOnm, with 3000 impurity 
configurations. 



(<5[ri dw) ] 2 ) 1/2 and (SIT^] 2 ) 1 / 2 with the correspond- 
ing predictions of Eq. 13411 . The two agree reasonably 
well and, significantly, the approximate values repro- 
duce the decreasing behaviour of (d^Ta^] 2 ) 1 / 2 for 1/1+ 
large. The simple hypotheses contained in Eqs. (|32I - I33|) 
thus explain quantitatively the two essential features of 
(i5[X'j tlw ' > ] 2 ) 1 / 2 , which are an overall reduction with re- 
spect to (<5[X'i um ' ) ] 2 ) 1 / 2 and a loss of universality repre- 
sented by a slow decrease with increasing disorder. 



D. Non-monotonic A-dependence of incoherent 
relative magnetoconductance 

We now consider the relative intrinsic domain wall 
magnetoconductance, defined as 



(Ag/g) dv 



9u 



9dv 



fju 



(35) 



For disordered systems this is a more useful quantity than 
the magnetoconductance Ag discussed in Section If 11 Bl 
as it measures the relative importance of the magneto- 
conductance effect inside the wall. In this section we con- 
sider only the incoherent regime; for this case, (Ag/g)d w 
is independent of L y since (Agi nc ) scales linearly with 
transverse system size. 

In Figure |Slwe plot ((Ag/g)d w ) as a function of 1/1+ 
(in units of 1/A). Various domain wall widths (p-p = 1,5) 
and up/down scattering ratios (p = 2,4) are shown. We 
see that ((Ag/g)^) is a monotonically increasing func- 
tion of disorder. This shows that the decreasing be- 
haviour of (A<7i nc ) for large disorder is due simply to 
the overall decrease in conductance with disorder. The 
monotonicity of ((Ag/ g) dw ) occurs because the amount 
of spin-mistracking, and hence also the enhancement of 
back-scattering, is proportional to the amount of impu- 
rity scattering, as discussed in Section Till Bl It can also 
be seen from Figure |H| that ((Ag/g)d w ) increases with p 
and decreases with pp, which is consistent with the ex- 
planation of the intrinsic magnetoconductance given in 
Section mm 

We note that the models of Refs. and assume 
diffusive transport, so that the results obtained therein 
correspond to the regime of large 1/1+ in Figure [S] The 
magnetoconductance predicted in these works depends 
on the up/down scattering ratio p, but is independent of 
1+ (provided it satisfies 1+ <C A). This is consistent with 
the levelling off of the curves in Figure [S] with increasing 
1/1+, and suggests that ((Ag/g)d w ) approaches a finite 
limit for large disorder. 

Figure El shows ((Ag / g)dw) as a function of wall width 
A, for p — 2, 4 and 1+ = lOnm, 40nm. For large A, we 
find a 1/A dependence in agreement with Refs. and 
HI . This can be understood because with increasing wall 
width the amount of spin-mistracking is reduced. In the 
limit A — > 0, on the other hand, the domain wall becomes 
an abrupt interface with no impurities and ((Ag/g)dw) 
becomes negligible, 51 ! In between these two extreme cases, 
there must therefore be a maximum value of ((Ag/ g) dw ) 
at a finite value of A. This value represents an optimum 
intermediate situation: a wall which is wide enough to 
contain an appreciable amount of impurity scattering, 
but narrow enough to cause significant spin-mistracking. 
From Figure [5] it can be inferred that this value of A 
decreases with increasing disorder (i.e. decreasing l a ). 
On the other hand, the different values of p affect the 
overall magnitude of ((Ag / g)dw) , but do not significantly 
change the position of its maximum value. 



14 



0.1 



J 

^ 0.01 



0.001 t 









p F = 1, p = 4 




Pf = 1, p = 2 - 








II 

!' .**' 

■i .*" 

t ■' 



2 4 6 

1/1+ [1/A] 



10 





\ f / 


- - 


X v X X 

x X X x 

X x X 

x >< LEFT " x 


x * x x 
x x * x 

»" DwV 


X x x * 

* x x 

x x x 

„ >< RIGHT ^ X x x 


x x x x 
X X x 

x X X X x 


X 

x X 
X » » « 


x v 
x X X * X x 

X X X X 




A ► 


' lad 



FIG. 10: Diagram illustrating the physical region considered 
in Section lTVI On either side of the domain wall are uniformly 
magnetized regions of length l s d representing the distance over 
which electrons propagate before spin relaxation. 



FIG. 8: Intrinsic relative domain wall magnetoconductance, 
{(Ag / g)dw) , as a function of disorder strength, 1/1+ (in units 
of 1/A), in the incoherent case. The different curves illustrate 
the dependence on wall width (pp = 1,5) and spin dependence 
of scattering (p = 2,4). The arrows indicate the evolution of 
the curves upon increasing these parameters. Other system 
parameters are as for Figure 0] 
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FIG. 9: Intrinsic relative domain wall magnetoconductance, 
{(Ag / g)dw) , as a function of wall width, A, in the incoherent 
case. Several values of mean free path (1+ = lOnm, 40nm) 
and spin dependence of scattering (p — 2, 4) are shown. The 
arrows show direction in which the curves evolve upon in- 
creasing these parameters. The two values of wall width cor- 
responding to Figure |H| (pf = 1,5) are indicated by dotted 
vertical lines. Other system parameters are as for Figure 0] 



IV. DOMAIN WALL IN A NANOWIRE: 
EFFECT OF SCATTERING IN THE LEADS 

To understand the effect of a domain wall on electron 
transport in a nanowire, it is necessary to consider not 
only the intrinsic effects described in the previous sec- 
tion, but also extrinsic effects arising from scattering in 
the regions adjacent to the wall. In a previous work 16 
we analysed such effects in the special case of a ballistic 
domain wall, using a circuit model which incorporates 
spin-dependent scattering in the regions on either side of 
the wall using a generalization of the two-resistor model 
of GMRiii The essence of this model is to treat transport 



of electrons with different spins in the adjacent regions 
as independent over a distance equal to the spin-diffusion 
length, Z sc j. This is the length scale over which the distri- 
bution between the two spin sub-bands is equilibriated 
through spin- flip scattering processes. In Ref. |n| trans- 
port in the adjacent regions was treated using classical 
spin-dependent resistors. Such an approach is valid un- 
der the assumption that the phase coherence length, L^, 
satisfies A < <C lsd- 

In this section we follow a similar method to incor- 
porate extrinsic scattering (outside the domain wall re- 
gion) into our model. However, in place of classical 
spin-dependent resistors, we treat the impurity scatter- 
ing in the regions adjacent to the domain wall using the 
same delta function scatterer model as inside the domain 
wall. The resulting model system is illustrated in Fig- 
ure This approach differs with respect to Ref. [l6| 
in that the adjacent regions are combined with the do- 
main wall using the scattering matrix combination for- 
mula (Eqs. @), rather than combining spin-dependent 
conductances using Kirchhoff's rules. We are interested 
in the limit l± <C / s d> for which the difference between 
these two approaches is not significant. However, the 
present approach has the advantage that it remains valid 
in principle for any value of l±/l s d- I n addition, phase 
coherence lengths > A can be considered within this 
approach (although we do not do so in this paper). 

As in Sect ion ITTT1 we consider conductance through the 
domain wall in both coherent and incoherent cases. In 
the coherent case we introduce a finite = A. The 
system is thus partitioned into phase-coherent segments 
of length which are combined incoherently. For the 
incoherent case, our approach assumes a vanishing phase 
coherence length, so no partitioning is necessary. 

The physical system in Figure EH has total length 
^wirc = 2^d + A. We will refer to this system as "the 
wire" since, by assumption, the transport outside this re- 
gion is in equilibrium between the two spin channels and 
therefore does not contribute to magnetoconductance ef- 
fects. Our main interest will be the relative magneto- 
conductance for the wire, (Ag/ <?) w ire, which is defined 
analogously to Eq. I|35|) by changing the intrinsic domain 
wall conductances to combined wire conductances. 
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Domain wall magnetoconductance in a wire: 
effect of disorder inside domain wall 



In Figure ITU fthick lines) we show ((Ag/g) wirc ) as a 
function of domain wall width A for both coherent and 
incoherent cases. To illustrate the dependence on dis- 
order, two values of mean free path are shown (Z + = 
lOnm, 40nm). In all cases we fix p = 2 and l s d = 
results for other values of these parameters are quali- 
tatively unchanged. For comparison with our previous 
work^t in which the domain wall was treated as ballis- 
tic, Figure^] a l so shows the magnetoconductance for the 
equivalent wire with no scatterers inside the domain wall 
(thin lines). 

In the narrow wall limit, A — > 0, the coherent regime 
is the physically relevant one. Figure ITTk shows that the 
coherent ((Ag / g) W ij- c ) (solid lines) tends to a maximum 
value as A — > 0. This is in contrast to the behaviour 
of the intrinsic magnetoconductance (Ag) studied in the 
previous section, which (for both coherent and incoher- 
ent cases) goes to zero in this limit. From our previ- 
ous workji& this result is expected since at A = the 
wall corresponds to a GMR interface^ for which there 
is complete mistracking of spin. It is interesting to note, 
however, that for small but non-zero A, ((Ag / g) w - 1IC ) is 
reduced for a coherent wall with disorder (Figure ITTk . 
thick lines) compared to a ballistic wall (Figure ^TJ thin 
lines) . Therefore, when scattering in the surrounding re- 
gions is taken into account, the presence of disorder inside 
the wall reduces the magnetoconductance effect of narrow 
domain walls. 

In the incoherent regime, the magnetoconductance is 
significantly reduced with respect to the coherent regime 
in the limit A — > 0. Furthermore, it is a non-monotonic 
function of disorder in this limit, having a maximum 
value at a wall width A > 0. In fact, this behaviour 
is an artifact arising from the lack of complete mistrack- 
ing when A — ► for transmission through an incoherent 
wall, which we discussed in Section III B II As we men- 
tioned above, in this limit the coherent regime describes 
the physical situation, so these effects are not physically 
relevant. 

In the wide wall limit A — > oo, the behaviour of 
((Ag / g) W i rc ) is qualitatively similar to that of ((Ag/g)^ w ) 
of Section III. The coherent magnetoconductance is sig- 
nificantly reduced with respect to the incoherent one and 
becomes negative for large A (in Figure ITTk this is most 
evident for l + = lOnm). These effects are due to the 
intrinsic properties of coherent transport through a dis- 
ordered domain wall which we discussed in Section [ill Bl 
For large A the incoherent regime is the physically rele- 
vant one. In this case the magnetoconductance decreases 
as 1/A 2 , which is the same parametric behaviour as the 

we will see in 



intrinsic case (Section IIII D|) . However, 
the next section that this effect is larger in magnitude 
for ((Ag/g) w - 1Ie ). In this case, comparison between the 
incoherent disordered wall (Figure ITTb . thick lines) with 
an incoherent ballistic wall (Tigure lTTb . thin lines) shows 
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FIG. 11: Magnetoconductance of a domain wall in a 
nanowire, ((Ag/ g) w i IC ) , as a function of wall width A for 
(a) coherent and (b) incoherent transport through the do- 
main wall region. Two values of mean free path are shown, 
1+ = lOnm, 40nm (thick lines). The magnetoconductance of 
an equivalent disordered wire with a ballistic domain wall (i.e. 
no impurities inside the wall) is also shown (thin lines). The 
inset in the upper figure shows the coherent ((Ag / <?) W ire) in 
the small A region. Other parameter values are L y = 5nm, 
p = 2, and i s( j = 4Z+, with iVs = 3000 impurity configurations. 



that the presence of disorder inside the wall enhances the 
magnetoconductance effect in wide domain walls. This 
enhancement has two causes. Firstly, for a disordered 
wall there is an intrinsic contribution due to scatter- 
ing from impurities inside the wall, which is absent in 
the case of a ballistic wall. Secondly, the disordered 
wall gives rise to an enhanced spin-mistracking (Section 
IIII A|) . which leads to an increase in the GMR scattering 
in the surrounding regions. 

To our knowledge no experimental data exists compar- 
ing the conductance through domain walls in wires with 
different amounts of disorder; the findings of this section 
suggest that such a difference may be observable. 



B. Enhancement of domain wall 
magnetoconductance from extrinsic scattering 

The total magnetoconductance of a wire containing a 
domain wall, (Ag/ g) wn . ci arises from scattering effects 
which may be classified as either intrinsic or extrin- 
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sic. Intrinsic effects are those due to scattering occur- 
ring within the region of the domain wall itself, and 
are characterized by the intrinsic magnetoconductance 
(Ag/g)d w studied in Section Till Dl Extrinsic effects, on 
the other hand, occur in the regions adjacent to the do- 
main wall, and arise from a GMR-like enhancement of 
back-scattering due to the spin-mistracking of electrons 
scattered by the wall. 

We now determine the relative importance of intrin- 
sic and extrinsic effects in the magnetoconductance of 
a domain wall in a nanowire. In order to do this, we 
compare (Ag/ g) w i TC , which contains both intrinsic and 
extrinsic effects, with (Ag/g)d w , which contains only in- 
trinsic ones. 

In the narrow wall limit, A — » 0, the amount of im- 
purity scattering inside the wall goes to zero, and hence 
((Ag/g)d w ) becomes negligible. On the other hand, Fig- 
ure shows that, in the physically relevant case of co- 
herent walls, ((A<?/g) w i, c ) achieves its maximum value 
in this limit. It is therefore clear that for narrow walls, 
extrinsic scattering effects are the dominant mechanism 
underlying the domain wall magnetoconductance. 

For general A, however, the situation is more com- 
plicated. We saw in Section IlII Dl that ((Ag/g)d w ) is 
significant for disordered walls. At the same time, the 
scattering from impurities inside the wall enhances the 
spin-mistracking of transmitted and reflected electrons 
(Section IIII Afl . which will lead to extrinsic scattering. 
We therefore expect that both intrinsic and intrinsic ef- 
fects will be significant in general. To determine the rel- 
ative magnitude of the two, it is necessary to compare 
((Ag/g) wirc ) and ((Ag/g) dw ) quantitatively. 

For this comparison we need to "renormalize" 
(Ag/g) dw to describe the same total system size as the 
wire, i.e. L W irc = 2Z s d + A. To do this, we calculate sepa- 
rately the total conductances (summed over all spin and 
transverse channels) for each of the three regions of the 
wire: the domain wall (<?dw), and the two surrounding 
regions (gioft and (/right)- Taking the incoherent com- 
bination of the three regions^ we then obtain a total 
conductance <7;„t over the length of the wire: 



1 - ffii 



1 — <?left . 1 — 5dv 



1 



Sright 



Slcft 



fjdv 



j;lit 



(36) 



Combining conductances in this way is equivalent to 
making the assumption that the two spin channels are 
equilibriated at the wall boundaries rather than at a dis- 
tance Z s d from the wall. 

By comparing g- mt with the conductance of a uniform 
region of length 2l s d + A, we calculate an "intrinsic" mag- 
netoconductance (Ag/g)- mt . This quantity represents the 
desired renormalization of (Ag/g)d w to the length L W1IC . 
We note that for ! + <A this can be calculated approxi- 
mately as 



(As/flOint = 



A 



2/ s d + A 



(Aff/ff) c 



(37) 
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FIG. 12: Ratio of incoherent wire and intrinsic domain wall 
magnetoconductance {(Ag / <?) w irc) / ((Ag / g)i n t) as a function 
of A, for illustrative values of 1+ and l s d- The system param- 



eters are p — 2 and L v 



lOnm. 



In Figure IT21 we show the ratio of ((Ag/g) w i IC ) and 



((Ag/g)i n t) as a function of A for several values of p 
and Z s d, in the incoherent case. We see that the ratio 
is largest in the limit of small A. As we mentioned pre- 
viously, this is because ((Ag/g)d w ) becomes negligible in 
the limit A — > while ((Ag / g) wilc ) becomes large. We 
note that the ratios shown in Figure lT2l apply to the in- 
coherent regime. As we saw in Figure ^2 ((Ag/g) w i re ) 
is much larger near A — > for the coherent case than for 
the incoherent one. Therefore, in the coherent case the 
corresponding enhancement of ((Ag / #) W ire) with respect 
to ((Ag/g)i nt ) will also be larger. 

For large A, Figure ^| shows the less obvious result 
that ((A<7/<7) w i re ) is substantially larger than ((Ag / g)- mt ) 
even in the case of wide walls (A ~ 200nm). The ratio 
depends on the parameter values, and for those of Figure 
1121 it is between ~1.2 (small l sd ) and ~2.2 (large l sd ) at 
A = 200nm. We thus see that extrinsic spin- dependent 
scattering effects are quantitatively important for domain 
wall magnetoconductance, even for wide walls. This is an 
important point, since the most widely accepted existing 
models^i&Si consider conductance only through the re- 
gion of the wall, and hence ignore extrinsic effects. 



V. CONCLUSION 

In this work we have studied spin-dependent electron 
transport through domain walls in disordered ferromag- 
netic nanowires. Using a novel numerical scattering ma- 
trix formalism, we considered a large range of impurity 
concentrations, corresponding to transport regimes from 
ballistic through to diffusive. In addition, our approach 
allowed us to consider spin-dependent impurity scatter- 
ing, and to distinguish between the fully coherent and 
fully incoherent transport regimes. 

The contribution of a domain wall to the resistance of 
a wire is due to two main mechanisms. The first one 
is the reflection of electrons from the domain wall itself 
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and impurities within the domain wall region. The sec- 
ond is related to the spin-dependent resistivity of the wire 
regions surrounding the domain wall. We therefore con- 
sidered first the spin-dependent "intrinsic" transport co- 
efficients of the disordered domain wall region itself, and 
then showed in the last section what can be expected for 
the resistance of the whole wire, which is the experimen- 
tally measurable quantity. 

The intrinsic effects stem from the presence of impu- 
rities inside a domain wall, which causes electrons to 
be scattered both with and without spin reversal. This 
leads to an enhancement of backscattering caused by ex- 
posure of electrons in a superposition of spin states to 
spin-dependent scattering. In the incoherent regime this 
gives rise to a reduction of the intrinsic domain wall con- 
ductance. The maximum reduction occurs at a disorder- 
dependent value of the domain wall length. 

In the phase-coherent regime, the relative amount of 
reflection with spin reversal is reduced with respect to 
the incoherent regime. For large disorder, this acts to 
suppress the backscattering enhancement observed in the 
incoherent case. Indeed, the overall effect in this regime 
is an enhancement of transmission through the domain 
wall. This can be understood as arising from a suppres- 
sion of weak localization by the domain wall, which is 
consistent with the findings of Ref. [^(J However, we 
point out several reasons why such a behaviour would 
be difficult to observe experimentally. 

The mesoscopic character of the transport in the coher- 
ent regime is underlined by a study of the conductance 
fluctuations. For a domain wall it is found that these are 
reduced with respect to the uniformly magnetized case, 
as a result of decorrelation between the components of 
a given incoming spin state which are transmitted with 
and without reversal of spin. Furthermore, the sample- 
dependent fluctuations of the domain wall contribution 
to the resistance can be significantly larger than its av- 
erage value in the diffusive limit. This should be kept in 
mind when experiments on individual samples approach 
the coherent regime. 

To address the experimentally measurable conduc- 
tance of a disordered nanowire containing a domain wall, 
we included in our treatment uniformly magnetized re- 
gions of length l s d (the spin diffusion length) on either 
side of the domain wall. Since the phase coherence length 
does not exceed Z sc j, the rest of the wire is taken into 
account as a classical resistance in series that does not 
contribute to the magnetoconductance. 

The spin-dependent scattering in the regions adja- 
cent to the wall results in an enhancement of the re- 
sistance. This arises from processes in which the spin 
of the electrons does not follow the domain wall mag- 
netization while they propagate through the wall. The 
spin-dependent scattering in the surrounding regions 
contributes an "extrinsic" domain wall resistance effect 
which is similar in nature to giant magnetoresistance. 
This mechanism is distinct from the intrinsic effect oc- 
curring inside the wall. 



As compared to the intrinsic domain wall resistance 
alone, the total resistance due to a domain wall in a 
nanowire is significantly enhanced by the extrinsic scat- 
tering. The enhancement is most dramatic for narrow 
walls, where the intrinsic resistance becomes negligible. 
However, it remains important (a factor of 2 for typical 
parameters) even for wide walls. 

In our previous study based on ballistic domain walls in 
a circuit modeU£ we found results in order-of-magnitudc 
agreement with a recent experiment^ The relevant 
regime for this experiment is intermediate between the 
coherent and incoherent regimes, and corresponds to do- 
main walls far from the adiabatic limit of zero spin- 
mistracking. In the present work, we have found that the 
presence of disorder inside the wall in this regime leads to 
a quantitative reduction of the magnetoconductance ef- 
fect as compared to the ballistic wall used in Ref. [l6L but 
is still within an order of magnitude of the experiment. 
This is acceptable given the uncertainties over various 
physical parameters relevant for the experiment (such as 
the precise value of the elastic mean free path). On the 
other hand, in the case of wide walls in the incoherent 
transport regime, we have found that disorder enhances 
the magnetoconductance effect. This regime is relevant 
for experiments in materials such as nickel and iron. Our 
model should also be applicable to the understanding of 
recent measurements of domain wall resistance in ferro- 
magnetic semiconductors. 43 

In both narrow and wide wall regimes, our results show 
that disorder inside and outside a domain wall leads to 
important effects in the magnetoconductance. Since dis- 
order is unavoidable in mesoscopic wires, our conclusions 
should be testable experimentally. 
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APPENDIX A: SCATTERING MATRIX OF 
DOMAIN WALL INTERFACES 

In this appendix we present formulas for the calcula- 
tion of the scattering matrices of the domain wall inter- 
faces s' L ' R ). These are obtained by forming scattering 
state solutions from the basis states of Eqs. |J5J and i|17|) 
for each side of the interface and matching the functions 
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and their derivatives at z — (left interface) or z — A 
(right interface). For the domain wall interfaces the ma- 
trix elements corresponding to scattering between differ- 
ent transverse modes, n' ^ n, are zero since the domain 
wall profile depends only on the longitudinal coordinate. 



The amplitudes f V±^. an and r^ an , 

satisfy 4x4 sets of equations. The equations for f^.^ 

and i^n-o-n rea d (for clarity we omit the n and L labels) 
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while for and i'^> we have 
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We present only the equations for the elements of 
explicitly; because of the symmetry of the spiral domain 
wall profile (Eq. H14I) ). the elements of can be found 
from Eqs. (|A1IIA2(1 after interchanging r a i a r' a , a , 
tg'g t'g'g and replacing k x -> -k\. 



APPENDIX B: DELTA FUNCTION 
SCATTERING MATRICES IN UNIFORM AND 
DOMAIN WALL POTENTIALS 

We now consider the scattering matrices for a delta 
function scatterer located in regions of uniform and ro- 
tating magnetization, s^(y a ) and s^ s \y a ). In the uni- 
form case, s^(y a ) can be obtained using the result for 



a spin-independent delta function potential^ The indi- 
vidual amplitudes are given by 



AS) 

cr'n' -.an 



Cl;gn = <W [Mg]^ n , (Bla) 

r {S) = r l{S) = t (S) - 5 , (Bib) 

g'n f ;gn g' n' \gn g'n'\gn n m \ ) 

where j\Ag is an N a x N a matrix with elements 

[Mg] n , n = S n , n + ■ (B2) 

/6 Y L! gn / u gn 

The scattering matrix for a delta function located in- 
side a domain wall, s^ s \y a ), is more complicated because 
the domain wall basis states ipg n (?) couple up and down 
components. Defining (3 n = A an A^g n we have 
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where T a ' a and are matrices of dimension N a > x N a and iV CT x _/V CT respectively, defined by 
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The amplitudes in Eqs. IjBll) and (|B3|I determine the matrices s^\y a ) and s^ s \y a ) that are composed numerically 
for each impurity configuration in order to obtain the total transmission coefficients in various regimes. 



APPENDIX C: DELTA FUNCTION SCATTERING IN WEAK DISORDER LIMIT 



In the limit of weak disorder, u a — > 0, transport through the impurity potential V(r) (Eqs. J22H2J) can be treated 
within the Born approximation. This approach yields the scattering amplitudes to lowest order in u a . For the 
uniformly magnetized case, we have 
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In the case of a domain wall, the situation is more com- 
plicated because scattering from the interfaces at z — 0, A 
leads to mixing of up and down spin channels. It is 
most convenient to first calculate amplitudes for trans- 
port through the disordered region < z < A, not includ- 
ing the interfaces, which we write £^*f„ and f^**^. 



Formally, these coefficients correspond to a scenario in 
which the rotating potential 9(z) — ttz/X is valid for all 
z, so that the asymptotic states are the domain wall ba- 
sis states ipan(r)- Within the Born approximation, these 
amplitudes are given by 
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To calculate the average scattering probabilities, it is 
necessary to take the squared magnitude of the quantities 
in Eqs. (|C1HC2|) and then average over the impurity posi- 
tions f a . Because the positions of different impurities are 
uncorrelated, the average of terms involving two different 
impurities are a factor l/kpL z smaller than those involv- 
ing a single impurity. They are therefore negligible, and 
the scattering is dominated by the single impurity scat- 
tering events^ 

We note that the probabilities corresponding to the 

amplitudes talari & n d t£n\an \ which are equal to unity 
in the limit u a — 0, cannot be calculated by simply taking 
the absolute square of the amplitudes in Eqs. (|Clh ') and 
(|C2b ). This is due to presence of the factor S n i n , which 
means that the term of 0{u^) in the amplitude (which 
is not given in Eqs. (|C1 h ) and l|C2b '0 also contributes 
a term of 0{u a ) in the corresponding probability. It is 
therefore necessary either to calculate the amplitudes to 
0(u^.) or, more conveniently, to use conservation of prob- 
ability to express it in terms of the other probabilities. 

With these considerations, we find the average proba- 
bilities in the uniform case as follows: 
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Summing over all n and n' , we can express the total 
spin-dependent transmission and reflection in terms of 
the mean free path l a (Eq. 
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where 
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Here Fi(N a ) is a form factor which takes account of the 
finite number of channels and reduces to unity in the 
limit Na — * oo. 



In the domain wall case, the probabilities correspond- 
ing to the amplitudes in Eqs. (|C2|) are given by 



an' u an 
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In general, it is non-trivial to incorporate the interfaces and thus obtain the total transmission probabilities for the 
domain wall, since if there is non-zero reflection from the interfaces then multiple paths through the disordered region 
need to be considered. The situation simplifies somewhat in the wide wall limit of Eqs. (|20|l . since to 0(l/p|) the 
interface reflection probabilities are zero. The transmission through the entire domain wall, including both interfaces 
and the disordered region, can then be found by multiplying the transmission of each part. Furthermore, since the 
interface transmission is diagonal in transverse channel number, only intermediate spin states need to be summed 
over. 

It is simplest to consider the incoherent case, since the probabilities are taken before combining with the interfaces, 
meaning that the results of Eqs. (|C3HC6|I can be used directly. In particular, for the transmission and reflection 
probabilities from the left we have 
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The total spin-dependent probabilities, (T^f^ ) and (R^™ ), may be found by summing over the previous ex- 
pressions, yielding quite cumbersome expressions. Of particular interest, however, is the transmission with spin- 



mistracking, (T^J), which is given by 
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Here F2(N a ) and F3(N <T ) are form factors which reduce to unity in the limit L y — > oo. Note that in the approximation 

used here we have fc + .p = fc-.F = &f and N+ = Eq. (|C8|I shows that the initial slope of T_™} as a function of 
1/1+ is always negative, which explains the observed behaviour the inset of Figure 0] 

Another important insight permitted by the perturbative approach is the spin-dependent difference in transmission 
AT± . Here we find 



(AT a ) 



71" V\ 



kpLy 



-F 3 (N a ) 



2N„ 
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Eq. JUTTll shows that (AT a ) increases linearly with 1/1+ for weak disorder, with a slope whose sign depends on a. 
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This explains the behaviour of (AT a ) observed in Figure 

Interestingly, the sum (AT+) + (ATI) in Eq. (RTiTI 
is zero, implying that the intrinsic magnetoconductance 
(Ag) is zero to first order in 1/1+. This is in agree- 
ment with our numerical results for (Ag), which show 
a quadratic dependence on 1/1+ for small disorder. 

Finally, we point out that for spin-independent disor- 
der, p = 1, Eq. jCTT)) predicts that (AT CT ) = 0. This is 
because in this case, the only spin-dependence of scat- 
tering comes from the difference in the up and down 
wavevectors, which is treated as zero in the approxima- 
tion of Eqs. Jl)|). 



APPENDIX D: SPIN-MIXING IN A 
ONE-DIMENSIONAL MODEL 

The results of Section UlI Al for the transmission and re- 
flection through a sequence of impurity scatterings show 
that repeated scattering from delta functions in a domain 
wall leads to a relative increase of the spin-mixing trans- 
mission and reflection coefficients, T- aa and ¥L- a<y . We 
now gain some insight into this mechanism through a sim- 
ple toy model which captures some essential features of 
the problem. In particular, consider a one-dimensional 
sequence of idealized spin-mixing scatterers for which 
there is zero reflection and for which the spin-dependent 
transmission probabilities are 

T = T'= ^; + £+ x e ~ e _y (Dl) 

Here e± represent the off-diagonal scattering. Techni- 
cally, this model is actually an example of a Markov 



chain, widely studied in probability theory^ 

The absence of reflection allows us to calculate the to- 
tal transmission through N[ scatterers by taking T 1 . Di- 
agonalizing T and taking powers of the diagonal compo- 
nents, we evaluate this as 

e+ + e- V e +l 1- X ) e + + e -X J 

where x = 1 — e + — e_. In the limit of infinitely many 
scatterers this reduces to 

T°° = ( 6 - e ~\ . (D3) 

We thus observe that the spin-flip scattering, when ac- 
cumulated over many events, leads to a complete mix- 
ing of spin, with the relative weightings determined by 
e CT /(e + + e_). In the special case of spin-independent 
scattering, e + = e_, all matrix elements in the r.h.s. of 
Eq. i|D3(l are equal and hence, for a given incoming spin 
direction, the transmission is mixed equally between both 
spin directions. 

The transmission T°° reproduces qualitatively the be- 
haviour of the fractional transmission from Section flll Al 
i.e. T a i a / (T aa + T-erer). Obviously, such a simplified pic- 
ture does not reproduce the actual results on a quanti- 
tative level. However, the qualitative resemblance sug- 
gests that it captures an essential feature of transmission 
through a disordered domain wall, namely, a spin-mixing 
effect which accumulates (and saturates) over successive 
scattering events. 
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